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ACOUSTIC THEORY OF AXISYMMETRIC MULTISECTIONED DUCTS 


By William E. Zorumski 
Langley Research Center 

SUMMARY 

Equations are developed for the acoustic field in a duct system which is made up of 
a number of connected circular and annular ducts. These equations are suitable for find- 
ing the acoustic field inside of and radiated from an aircraft turbofan engine. Acoustic 
modes are used as generalized coordinates in order to develop a set of matrix equations 
for the acoustic field. Equations for these modes are given for circular and annular ducts 
with uniform flow. Modal source equations are derived for point acoustic sources. Gen- 
eral equations for mode transmission and reflection are developed and detailed equations 
are derived for ducts with multiple sections of acoustic treatment and for ducts with cir- 
cumferential splitter rings. The general theory is applied to the special case of a uniform 
area circular duct with multisection liners and it is shown that the mode reflection effects 
are proportional to differences of the acoustic admittances of adjacent liners. A numeri- 
cal example is given which shows that multisection liners may provide greater noise sup- 
pression than uniform liners. 


INTRODUCTION 

Intense noise is radiated from modern turbofan aircraft engines. This noise may be 
broadly classified as externally generated noise and internally generated noise. The main 
external noise source is the engine jet. The main internal noise source is the turboma- 
chinery stages. Steady and unsteady flow through these stages causes fluctuating loads on 
the blades that generate an acoustic field in the engine. This acoustic field is propagated 
away from the blades, transmitted through the engine ducts, and radiated to the far field. 

Research in the area of internal aircraft engine noise has been traditionally divided 
into three fields: turbomachinery noise generation, sound transmission in ducts, and 
radiation from ducts. Although extensive literature is available in each of these fields, 
little published work which integrates these fields into a general theory exists. A first 
step toward organizing these fields was taken by Lansing and Zorumski (ref. 1) who treated 
the problem of a uniform multisectioned duct with a matrix approach. Alfredson (ref. 2) 
has used a similar approach to solve the problem of a hard-walled acoustic horn by match- 
ing the sound fields in a number of disk-shaped elements which approximate the horn shape. 



Matrix concepts have been used advantageously in other engineering disciplines, such as 
structures and electrical networks, to integrate and organize complex problem areas. 

The electrical theory has already been carried over into acoustics in the one-dimensional 
approximations used in muffler theory (ref. 3); however, no accepted method has appeared 
previously for systematically treating multidimensional problems in duct acoustics. 

This paper develops a detailed theory for interactions of the acoustic fields of a 
number of interconnected simple geometric regions. All acoustics problems involving 
sound generation in, transmission through, and radiation from duct systems are shown to 
be reducible to the solution of a matrix equation [wJ{a} = { Q }. The matrix formulation 
is suitable for numerical solution on modern digital computers. The matrix notation also 
gives physical meaning to the numerical operations used in determining the acoustic field. 
For convenience, a number of simplifications are made in modeling the real duct systems; 
however, these restrictive simplifications may be removed without altering the essential 
structure of the resulting equations. The acoustic field is assumed to be governed by the 
linear convected wave equation. The geometry of the ducts is approximated by a number 
of interconnected circular and annular cylinders. The steady flow field is approximated 
by a uniform velocity along the duct axis with no temperature or density gradients. The 
acoustic treatment on the walls of the ducts is assumed to be acoustically linear. Point 
sources of sound are used in place of actual distributed sources, and the equations which 
couple the internal field to the external radiated field are assumed to be known. The 
matrix formulation is capable of representing more realistic physical situations than the 
one described above and will be useful as a standard for integrating various theories into 
a common framework. 

Another important purpose of this paper is to show that the interactive effects of 
acoustic treatment in different regions of an aircraft engine may be used to reduce the 
total noise radiated from the engine. Past efforts to reduce aircraft noise by using 
acoustic treatment have depended on the attenuation of waves in ducts. An equally impor- 
tant consideration, wave reflections at impedance discontinuities, has not been utilized. 

It was shown in reference 1, however, that wave reflection effects are important and that 
these effects may be used in suppressing noise transmitted through ducts. The matrix 
formulation adopted in this paper allows these effects to be easily incorporated in noise 
reduction studies. 


SYMBOLS 


A wave amplitude 

a = fea* 
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inner and outer radii, respectively, of an annulus 


a*,b* 
b = kb* 

c a 

D/Dt 

D/Dt* 

D(S2) 

W), 1 

D(sin (f>)J 
E(Ar) 

f 

1 

T* 

r 

J,Y 

K™(b,S2) 

b 

M 

N 

P 


ambient speed of sound 

dimensionless comoving derivative based on steady flow convection (eq. (6a)) 
comoving derivative 
eigenvalue determinant 

radiation directivity factors 

eigenfunction for annular duct (eq. (30)) 
unit vector in z-direction 
frequency, Hz 

scaled complex body force amplitude (eq. (4d)) 
body force 

eigenfunction product integral 

Bessel functions of first and second kind, respectively 

infinite integral (eq. (64)) 
wave number, cu/c 
uniform flow Mach number 

eigenfunction normalizing factor (eqs. (33) and (34)) 
scaled complex acoustic pressure amplitude (eq. (4b)) 
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P m (r,n) 

p* 

Q 


R 

R = |r| 
R = bR* 
R* 
r,8,z 

T 


t = U>t* 


t* 

V* 

v 

W 

0 

7 

A r 

6(x) 


axial Fourier transform of circumferential harmonic of p 
pressure 

volumetric source strength, dimensionless, or source mode amplitude when 
subscripts and superscripts are used 

reflection coefficient, used with subscripts and superscripts 


position vector 
cylindrical coordinates 

transmission coefficient, used with subscripts and superscripts 

time 

velocity 

scaled complex acoustic velocity amplitude (eq. (4c)) 

duct wave coefficient, used with subscripts and superscripts 

acoustic admittance, dimensionless 

ratio of specific heats of gas 

boundary condition operator (eq. (24)) 

Dirac delta function 
Kronecker delta function 
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€ « 1 


small scaling parameter 


A dimensionless radial wave number 

p scaled complex acoustic density amplitude (eq. (4a)) 

p„ ambient density 

Or 

p* density 

i//(r) normalized eigenfunctions for annular duct (eqs. (31) and (32)) 

ft axial Fourier transform parameter (propagation constant) 

ca> circular frequency, 27rf, radians/sec 

V = V*/k 

V* divergence, or "del" operator 

9/9n partial derivative normal to boundary 

Subscripts: 

a=fea* 

b=feb* 

m circumferential mode number 

o denotes source 

p. radial mode number and indicates row number in matrix 

v radial mode number and indicates column number in matrix 

Superscripts: 

j,k indicates plane number 
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l,k indicates general values of j and k 

+ upstream waves 

downstream waves 
Matrix notation: 

row matrix 
column matrix 
square matrix 

diagonal matrix 
identity matrix 


LJ 

o 

[] 

[] 

M 


GENERAL THEORY AND CONCEPTS 

This section will present the basic theory of approximating the acoustic field of a 
typical aircraft turbofan engine such as in figure 1. In order to bring the problem at hand 
within the scope of present analytical capability, the engine shown in figure 1 will be 
modeled by the hypothetical engine shown in figure 2. This engine model is simple 
enough for analytical purposes yet it retains the gross geometrical and flow character- 
istics of the real engine. The curved engine surfaces have been replaced by uniform 
cylindrical surfaces to permit a simple representation in cylindrical coordinates, the flow 



Figure 1.- Cross section of typical aircraft turbofan engine. 
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Figure 2 ,- Engine acoustic model. 

field has been replaced by a uniform flow with Mach number M, and the finite size turbo- 
machinery blades have been replaced by blades which have infinitesimal cross sections. 

Governing Equations and Boundary Conditions 

It will be assumed that all variations in the fluid state are isentropic, except in 
some small regions which may be considered separately, so that the energy equation may 
be given as 


Dp* _ TVl Dp* ^ 

Dt* p* Dt* 

In addition to equation (1), the continuity and momentum equations are required to deter- 
mine the fluid motion. These equations are 

2£L + p*v* • v* = o ( 2 ) 

Dt 


and 


p* 2X1. + V*p* = p*f* 
Dt* 


( 3 ) 


In theoretical acoustics, it is convenient to use dimensionless variables where the 

length and time scales are 1/fc and 1 / oj, respectively. Dimensionless lengths are 

formed by multiplying actual lengths by the wave number, and a dimensionless time is 

introduced by multiplying by the circular frequency u). The forms of the continuity and 

momentum equations may be preserved then if the density and pressure scales are p„ 

a 

o 

and p a c a and the velocity is referred to c a where p a and c a are the ambient den- 
sity and speed of sound, respectively. 
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The flow field is separated into a steady part, or mean flow, and a small amplitude 
periodic variation of frequency w. If a small scaling parameter € « 1 is introduced 
to account for the small amplitude of these variations, then the dependent field variables 
may be expressed in complex form as 


P* = p a ( 1 + ep) e _lt 

(4a) 

P* = P a c a^y + € p) e_lt 

(4b) 

V* = c a (-Me z + evj e _lt 

(4c) 

f* = cuc„ef e ^ 
a 

(4d) 


and the independent variables are 

t* = — (5a) 

(jO 

and 

r* = 1 (5b) 

fc 

Equations (4) show that the mean flow field is assumed to have constant density, 
pressure, and velocity. Note that equation (4c) implies that the mean flow is conven- 
tionally in the minus z -direction. 

With the approximations of equations (4) and the definitions of equations (5), the 
differential operator for the material derivative and "del" become 


— = J— + ev • V 
Dt* \Dt 


and 


v* = tv 


(6a) 


(6b) 
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In equation (6a), the part of the material derivative due to the acoustic velocity is 
shown explicitly so that it is understood that the dimensionless material derivative con- 
tains the convection effect of steady flow only. 

Substituting equations (4) and (6) into equations (1), (2), and (3) and collecting coef- 
ficients of like powers of € show that the acoustic equations are 

(7a) 

Dt Dt 

— + V • v =0 (7b) 

Dt 

— + Vp = i (7c) 

Dt 


The wave equation is derived by eliminating the velocity from the continuity and 
momentum equations. When the steady flow velocity is a constant, as in equation (4c), 

the operators — and V ■ commute; therefore, the wave equation becomes 
Dt 


v 2 P = s!E + v.r 

D Z t 


where 


D_ 

Dt 



( 8 ) 


( 9 ) 


Eversman and Beckemeyer (ref. 4) have shown that an approximation to the wall boundary 
condition where there is a sudden change in the mean flow velocity is the so-called parti- 
cle displacement boundary condition 


ig E?e.|E = o 

V 8n 


do) 


This boundary condition will be used later in the present paper. 
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Concept of Reflection and Transmission 


In this section, a brief heuristic discussion of transmission and reflection concepts 
is given. These concepts are developed from the physics of one-dimensional wave prop- 
agation. The equation forms postulated here guide the later development of a mathemat- 
ical theory of transmission and reflection. 


One-dimensional waves .- Figure 3 shows sche- 
matically a duct section where there is a wave repre- 
sented by the symbol A, traveling forward through 
the plane j. Due to some properties of the duct sec- 
tion between planes j and k, the amplitude of this 
wave may be modified before it reaches plane k. 

This change is shown symbolically by the equation 


A +k _ T +k+j A +j 




Figure Mode transmission and 
reflection in simple duct 
section. 


(ii) 


The factor T + ^ + i in equation (11) is called the transmission coefficient for the section. 
There may also be a wave reflected from the section, as shown symbolically by the 
backward- moving wave A 1 in figure 3. This reflection effect may be represented by 
the equation 


A 1 = R"i + i A +j (12) 

where R~l + 1, which is called the reflection coefficient, depends on the properties of the 
duct section between planes j and k. 

_ 1 ^ 

Similar effects may be expected if there is a backward- moving wave A incident 
on the section. In this case, the transmission and reflection effects are represented by 


A-i = T‘i‘ k A- k 

(13) 

A+ k „ R +k-k A -k 

(14) 


If there are sound sources in a duct 
section, as shown in figure 4, there will be 
waves propagated away from the source in 
both directions. This wave generation is 
represented by source terms Q + and Q 
for the forward and backward waves, 
respectively. 
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Figure 4.- Mode generation by source in 
duct section. 



The general sound field in a duct section is a superposition of the fields due to 
incident waves from both directions and due to sources within the section. Consequently, 
two general equations may be written to describe the effects of a section on the acoustic 
wave at either end of the section. 


A +k = T +k+ U + j + R +kk A -k + Q +k 

(15) 

A“i = T -j- k A' k + R-j + U + j + Q-i 

(16) 


Equations (15) and (16) are linear algebraic equations which represent the super- 
position of the sound fields due to various effects. Equations of this form may always be 
derived since the governing partial differential equations for the sound field are linear. 

Two-dimensional waves .- There are many waves with characteristic wave-front 
shapes which may propagate in a channel. A general two-dimensional wave in a channel 
is formed by superimposing these characteristic waves. The amplitude of this general 
wave is then a vector formed from the amplitudes of the characteristic waves. There- 
fore, the wave generation, transmission, and reflection effects must be represented by 
matrix equations rather than the simple scalar forms in equations (15) and (16). These 
matrix equations are assumed to be 



Some detailed examples of these equations will be presented in the next section. 

SOURCE IN INFINITE UNIFORM DUCT 


The sound field due to a monopole source in an infinite annular duct may be found by 
solving equation (8), with V • f replaced by a Dirac delta function, subject to the bound- 
ary conditions in equation (10). In cylindrical coordinates, the wave equation (eq. (8)) is 


a2 f + ? + \ ^2 ■ 2iM F + (1 ' M2 )^ + P = T 1 Q6(r " r o) 6 ( 0 -0o) 6(z-z 0 ) 

8r 2 8r r 2 dd 2 dz 8z 2 


(19) 
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and the boundary conditions are 


9p 

— ± i/3_ 
8r r 


i + M — 

3z 


= 0 


f r=b 
or 
r=a 


+ for r = b\ 
- for r = ay 


( 20 ) 


The solution to this problem in the special case where /3 r = 0 has been given previously 
by Drischler (ref. 5). Following Drischler's approach, it is assumed that the solution to 
equations (19) and (20) may be found as a circumferential series of inverse Fourier 
transforms 




dO 


m=-°° 


Substituting equation (21) into equations (19) and (20) gives 

Q6(r-r 0 ) -i(m0 o +Qz o ) 


1 _d_ 

r dp m (r,n) 

r dr 

| 

dr 




P m (r,n) = 


where 


( 21 ) 


( 22 ) 


2\o2 


A = \1 + 2MQ - (1 - M z )fi 


and 


:[p m ( r ^)] = 


dp (r,fi) 2 

m ± i|8 r (l + Mft) 2 p m (r,n) 


dr 


r=b 

or 

r=a 


- for r = b 
+ for r = a 


(23) 


(24) 


P m (r,S2) = AJ m (Ar) + BY m (Ar) - - Q e 


-i(me 0+ J2z 0 ) j J m(^o) Y m (Xr) (r ^ r 0 ) 


J m (Xr) Y m (Ar 0 ) (r ^ r 0 ) 


(25) 
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The first two terms of equation (25) are homogeneous solutions of equation (22), and the 
last term has a slope discontinuity which is introduced by the delta function. The coef- 
ficients A and B in equation (25) must be chosen to satisfy the boundary conditions of 
equation (24). Therefore, 



where 


Dm(fi) = A £ 


'm 


m 


- A, 


m 


'm 


(27) 


Substituting equation (26) into equation (25) gives the expressions for the transforms of 
the pressure harmonics p m (r,ft) in the regions r = r Q and r = r 0 as follows: 




-i(mfl 0 +nz 0 ) 


2 D m(«) 


- 


- A a[ J m] T ^m(' vr ) ^bj^mj Y nl (Ar 0 ) 

A b[ Ym ] " ^bpmj Y m (\r)J jA a |Y m j Jm(^ r o) “ 


( r £ r o> 
(r £ r 0 ) 


(28) 


It may be shown that D m (ft) and the other expressions on the right side of equation (28) 
are analytic functions of ft and that the inverse transform (eq. (21)) may be evaluated by 
the theory of residues. If z > z 0 , a contour which circles the upper half of the ft-plane 
is used and, if z < z 0 , a contour circling the lower half-plane is used. Substituting equa- 
tion (28) into equation (21) gives 



The constants ft m ^ are the zeros of equation (27), and equation (23) defines A m ^. 
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In equations (29) and (30) the + sign is used for z > z Q and the - sign is used 
when z < z G . It will be convenient in the following analysis to use characteristic func- 
tions ^ m/i ( r ) which are proportional to but which are normalized such 

that 



Consequently, * s defined by 


i 1/ (r) 

r m/i' ’ 



(32) 


where 


n 2 = r l 

m/i J„ 


rE ( \ r] dr 
ml m n 1 


The integral (eq. (33)) may be evaluated (ref. 6, p. 135), in closed form 

ib 

,41 ^ 


N 2 = 
m/i 2 ' 


1 - 


m‘ 


X 2 r 2 '^ 
A m/i r 


1 + Mn mnl 
2 \ m/ij 


'■m/i 


E m \n M r 


(33) 


(34) 


The characteristic functions 4 / mfJL ( r ) are not always orthogonal. If * \ mV , then 

it may be shown that if 


r b 

^m/iy " j ri ^m ^m^ r ^ 


(35) 


then 


m/ii’ 


2M + + o m „) 


1 - 


2M ♦ M*(0 m(1 + S2 ml , 


ia 0 a 1; m„ (a l , *W (a) + lb V»,l l l 


nijLi 


(36) 
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If M = 0 or if both j3 a and /3jj are zero, then equation (36) shows that 
^m( X m p r ) is orthogonal to that is, l m(lU = 0 when p * v. Equations (31) 

to (36) are valid for both upstream (+) modes and downstream (-) modes so that (±) super- 
scripts may be placed in these equations when they are required. 

Tester (ref. 7) has shown that N m ^ may be zero in certain cases which corre- 
spond to values of the wall admittance where = 0- These cases are excluded 

here. 

The expression for the sound field due to a single source may now be written in 
terms of the characteristic functions as follows: 

OO OO 

p(M,z) = l e ira0 l Q* m(t <z) ^(r) (37) 

m=-°° p=l 


where 



The expression for the modal source strength (eqs. (38)) may be extended easily to 
represent a number of superimposed monopole sources or a source distribution 

Q(r 0 ) over some finite region. The method of extending this solution to obtain a dipole 
source representing the force on a rotating blade has been given by Drischler (ref. 5). 

UNIFORM MULTISECTIONED DUCT (FIG. 5) 


In practical problems such as the uniform 
multisectioned duct depicted in figure 5 there 
are a number of finite duct sections which are 
interconnected. Therefore, the infinite duct 
solution (eq. (37)) must be generalized to 
account for the effects of these finite sections 
and relationships must be developed to account 
for the acoustic coupling between sections by 
(1) defining the acoustic field in terms of wave 



1 2 N-2 N-l N 

Figure 5.- Uniform multisectioned ducts. 
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amplitudes at a number of planes in the duct, (2) developing transmission and reflection 
relationships to relate these wave amplitudes, and (3) solving a matrix equation for the 
wave amplitudes. 


Field at a Plane 

In a uniform duct section containing the plane z^ = Constant, the acoustic field may 
be expressed as a sum of the acoustic modes Since the duct sections are axi- 

symmetric, the discussion may be limited to the circumferential harmonics of the sound 
field p m (r,z) which is given by the following equation: 

00 

Pm (r ' z)= 1 

11 = 1 


m /! ' r mp 


(r) e 

~ TY1 I 1 x / 


itt + ^ (z-z^) 
m /Li V / 


if2”j (z-z-i) 

- A' J (r) e ’ 

m/i r mju' ' 


(39a) 


The coefficients A + ' 1 and A' J in equation (39a) are the' amplitudes of the acoustic 

m jji m 

modes at t ). The negative sign in front of Aj^ has been introduced so that the mode 
reflection coefficients will have a positive real part. In order to represent the sound field 
everywhere in a complicated system of ducts, a large number of planes zl are intro- 
duced. The fundamental problem of determining the acoustic field is then reduced to find- 
ing the wave amplitudes A"^ and A"^ at these planes. Thus, if there are N planes 
then there must be 2N sets of equations relating these coefficients. Each of these sets 
contains an infinite number of equations; however, it is necessary in practice to truncate 

these sets to a finite number which determines a finite number of the coefficients A^| 
in equation (39a). Matrix notation is used to write these finite sets of equations. With 
matrix notation, equation (39a) becomes 


Pm (r ’ z) 


(r) 


m/i 






i//”j (r) 

■ TVl f f ' 1 


m/i 







(39b) 


Transmission and Reflection Matrices 

The second step in solving the multisectioned duct problem is to find the transmis- 
sion and reflection equations which relate wave amplitudes at the specified duct planes. 
Two methods of finding these equations will be illustrated in this section: the method of 
solving the wave equation in the section between two planes and the method of matching the 
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solutions at two planes which are close together by use of the continuity and momentum 
conservation conditions. 

Uniform duct.- The transmission and reflection matrices for a uniform duct section 
may be obtained by inspection of equation (39b). Since there are no reflections in a uni- 
form duct, 


R ±j± k 


0 

nij \iv 


L J 


(40) 


and since the change of the mode amplitudes is given by the exponential factors in 
equation (39b), 


I rp±k±j ! 

I mjuij 



\iv 


jj 


(41) 


k i 

In equation (41) the positive sign is used for z > z J 

k . i 
z < z J . 


and the negative sign is used for 


This example shows how a solution to the wave equation may be used to find trans- 
mission and reflection equations. It may be shown that, if a Green's function for a region 
exists, then transmission and reflection equations of the form of equations (17) and (18) 
may be derived from this Green's function. 


Admittance discontinuity .- When there is a dis- 
continuity in the wall admittance of an annular duct 
section with uniform inner and outer radii as shown in 
figure 6, then the conditions of conservation of mass 
and axial momentum at the interface may be used to 
derive equations for the transmission and reflection 
matrices. These conditions are satisfied if the axial 
velocity and the pressure are continuous at the inter- 
face. The axial velocities due to the waves are related 
to the pressures through the axial modal admittance 


matrix 


which the momentum equation shows is 


ix , 

+ Mfl mi; ) . Consequently, for planes t) and 

' J 

z^, which are an infinitesimal distance e apart as 
shown in figure 6, 


H c h- 

j k 

Figure r. . - Duct with admittance 
discontinuity. 
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I/T 3 (r) 


'fl-J 

p mv 




^ +k (r) 


q +k 


+k 

L mi^ 


^ +] (r) 


ft 


+J 

ini' 


+J 

‘mi' 


ib~ k (r) 


,-k IJ A -k 

mv 1 1 mv 


(42) 



* 


+k 

mv 





(43) 


The terms in equations (42) and (43) have been arranged so that the left-hand sides con- 
tain waves traveling out of the region between planes z^ and z k while the right-hand 
sides contain terms traveling into the region. In order to obtain an equation like equa- 


tion (17), it is necessary to eliminate from equations (42) and (43). In the special 

case where M = 0, this elimination may be carried out by utilizing the orthogonal prop- 
erties of the modes , r / ^^( r )- If M = 0, then premultiplying equation (42) by 


r (r)? dr, premultiplying equation (43) by r 
m/i. j 


m/i 


( r ) • clr ’ integrating each 


equation from a to b, and subtracting the resulting matrix equations would eliminate 

A"Q. The elimination process will be started in the same manner in the general case 

where the modes are not orthogonal so that the resulting equations will reduce to the 
orthogonal mode case. Premultiplying equations (42) and (43) by the modal vectors and 
integrating give 



\ 


- 



T-j'j 

m/i v 

i3“ j 

j^mvj + 

j-j+k 
m/i v 

o+k 

p mv 

{ A mv} 


” 



” — 

* 


T-j+j 

m/i v 

/3 +j 

fa,} • 

j-j-k 

m/iv 

/3 -k 

{ A m k v} 


(44) 


T j”j 
m/i v 



+ 


r-j+k 

mjitv 



t-J+J 
mp v 




(45) 
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where 


r - 


r ^ 


tIK 

X m \iv 

4 ■■ 

J a 


\L K (r) 



_ 



( 46 ) 


is called the mode integral matrix. 


Eliminating 



from equations (44) and (45) gives 


w -j +k 

(a*) . 

w-j + j„ 

[a + 0 + 

w -i- k 

f A_k ) 

mju v 

1 mv J 

mpf 

V mvj 

m/if 

l mV ) 


(47) 


where 


- 


r 

-1 

r -| 



* 

r- ^ 

-1 

- 

W 1 * 


l LL 


K 

(3 k 


& L 

jtt 


jL K 

mju^ 


mptt' 


A m jiv 



H m/i 

m/if 


l m\iv 

- 


- 


- 

J 



- 




(48) 


Now, a premultiplication of equation (47) by 


i-l 


W 


-3+k 

mfiv 


j 


gives an equation of the form of 


equation (17) from which the transmission and reflection matrices may be identified as 



r 

-1 

-j 

rn+k+j _ 

w -j +k 


w'j + j 

m iiv 

m/it' 


m/in 


(49) 


R+k-d . 

w -j +k 

m/un 

mji v 


7 . 


W"j' k 

muj'i 

J 


(50) 


It may be shown, by eliminating the vector 



from equations (42) and (43), that 


W 


+k-j | 
m/ar; 



w + k -k 
m [iv 



+ 


m /u^ 



(51) 
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and comparing equation (51) with equation (18) shows that 


- - 


r i 

-i 

r -» 

rp” j'k 



w + n 


w +k ' k 

mjui' 


mjjLV 




( 52 ) 


- - 


r n 

-1 

r -i 

R' j+j 



w +k_ i 


w +k+ [ 

m/iv 






(53) 


When M = 0 


so that the modes are orthogonal, the expression for the matrix 


W 


IK 

mij.v 


may be simplified since 


[“ 

mpy 


is an identity matrix in this case. 


Carrying out the inte- 


gration in equation (46) gives its elements 


I 


IK 

m/ii' 




(1 - M 2 ) 





+ Mft 



2 



and in the special case when M = 0, equation (54) reduces to 
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Therefore, if M = 0, 


W 


L K 

m py 




(b) „(b) ♦ 


+ Sl K 

mp my 


\b L (a) (a) 

^mp' ' ^my' ' 


(56) 


Equations (50), (53), and (56) show that the reflection matrices are explicitly proportione d 
to the admittance changes in the duct since W“j“ k and are proportioned to the 

admittance changes. Note that and are not proportional to the admittance 

changes since when M = 0. It may be seen from equation (46) that, since 

X mp = X mp if M = 0 > X mMy = I mpy = 6 py ; therefore, e q uation ( 48 ) shows th at 

W’j + j = 26 and W +k-k = 26 fi +k 

mpy py my mpy py my 


Radiation Equations 

If equations like equations (17) and (18) are repeatedly written for a number of 
planes, the number of unknowns always exceeds the number of equations. In order to 
have a complete set of equations, it is necessary to specify radiation conditions at bound- 
ary planes such as the engine inlet or exhaust duct termination planes. These radiation 
conditions may be formulated as special cases of equations (17) and (18) where the trans- 
mission term is dropped (no incoming wave) and the source term is dropped (no sources 
are exterior to the regions under consideration). Thus, the general form of a radiation 
equation is 



R -t+t 

mpy 



(57) 


The problem of determining the radiation reflection matrices is very difficult; however, 
some results are available in special cases. 

Flanged duct .- In reference 8, the reflection matrix is derived for the case of an 
annular duct without flow and with an infinite rigid flange. The sequence of analysis in 
this case is to use the radiation directivity factors 
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"\ 


mp 


(?) = { r 


Jm (Tr) • J m ( Tr ) 


m/i 


m/i ^ m/i 


nr 


m/i 


(58) 


(without flow, X mM = X^ = X mM , and -O^ = to calculate, first, the 

generalized radiation impedances by the infinite integral 


mpi' 


pOO 

= -i\ t(t 2 - 
Jn 


1) 


- 1/2 


D (t) 

mji' ’ 


D m*< T > 


dT 


(59) 


and, second, the radiation reflection coefficients by the matrix equation 
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m [iv 


t lr 



- - 


- - 



I 

+ 

Z mju^ 


_ 

- J 


J 



I 


r "i 

* 


z 

ft , 


m [iv 

mv 


- 




(60) 


The far-field acoustic pressure for the flanged duct may be calculated from the acoustic 
mode amplitudes at the duct termination and the radiation directivity factors by the 
summation 


p(R 



sin <fi 


(61) 


where 


V = ft 
zm/i m/i 


m/i 


+ A 


m/i 


(62) 


Unflanged duct .- A more realistic radiation model for aircraft applications, the 
semi-infinite unflanged duct, has been solved for the case of a hard- walled circular duct 
without flow in reference 9. In this problem, the radiation reflection coefficients are 
given by 
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R _ 1 b x mu 1 + n m/i 1 + fi mi' F m L 0 \ y m L 0 ' 

•w- - 2 r~M X 2 a +n K * b,s H K + r 

Y m n K ' b A m ^ - m m/i mp ini' ' ' ' ' 

where 


K™(b,z) = exp] ■-- 
+ 2m 


l 06 e - 2 K m W, 2 -b 2 /l^l^-b ,i l|d, 


77 - bz 


.2 ,2 


The radiation directivity factors for the unflanged, hard-walled circular duct are 


H 1 + 0(°m M + cos *) Wm M J m< b sin *) sin * 

D mp^ - T— 1 -4 (65) 

2(1 + cos 0)(\2 . sin 2^j K ™(b,cos 0) 


and the far-field acoustic pressure is given in terms of the mode amplitudes which are 
incident on the duct termination and the directivity factors by the summation 


p(R , 0 ,<P) 


I ‘ 


■( I m l +1 


A D (d>) 
m(x m ii yY> 


Application to Circular Duct 

Theory .- As an example of the application of the general theory, a circular duct is 
used. Flow is not considered since it has been shown that the only essential effect of flow 
is to make the modes nonorthogonal. The duct is shown in figure 7 with the left end 
closed by a rigid plate and the right end terminating in an infinite flange. The flange is 
not essential to the solution of this problem since unflanged radiation reflection coeffi- 
cients and far-field pressure formulas are available in this case. The rigid plate has a 

small orifice in the center through which there is an oscillating flow. The plane of the 

1 9 

plate is designated as z , the plane of the right end of the first section by z , the left 

end of the second section by z^ (z^ = z%), and so on to the plane of the flange z®. 

The velocity distribution at z* may be represented by the product of a Dirac delta 
function and a source strength Q. 
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Radiated 

sound 


Figure Source in a multisection duct. 


V(r) = ^ Q6(r) 


(67) 


Since the source and the duct are axisymmetric, the acoustic field is axisymmetric and 
only the circumferential mode where m = 0 is considered. The velocity distribution 
(eq. (67)) must be matched by the acoustic field at z 1 , consequently, the momentum equa- 
tion (eq. (7c)), the acoustic field equation (eq. (39b)), and the modal orthogonal conditions 
(eq. (36)) are used to show that 



Equation (68) is of the form of equation (17) with the transmission matrix dropped and 
represents a boundary condition on the acoustic field. 

The transmission and reflection matrices for the uniform sections (eqs. (40) and 
and (41)), for the wall admittance changes (eqs. (49), (50), (51), and (52)), and for the radi- 
ation reflections (eq. (57)) are combined with the source equation (eq. (67)) to give a com- 
plete set of equations for the wave amplitudes in the duct as follows: 
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In equation (69), only nonnull matrices are shown. 



The matrix equation (69) is equivalent to equation (23) of reference 1. In order to 
reduce equation (69) to the form of equation (23) of reference 1, it is necessary to elimi- 


nate the vectors 



from equation (69). 


This 


operation is easy since these vectors are not reflected in equation (69); that is, they are 
not multiplied by reflection matrices. 


Results.- This problem was considered previously in reference 1 and it was shown, 
by using idealized wall admittance values, that a duct with three liner sections with dif- 
ferent wall admittances can be made to radiate less acoustic energy than a duct with an 
optimum uniform liner. The wall admittances were chosen to increase the effect of the 
reflection matrices by having large admittance changes between adjacent sections. 


Based on these previous results, a set of duct liners were designed which would 
have large admittance changes at a specified frequency. The liners were made with single 
layers having different backing depths. The first and third liners had backing depths of 

3 1 

about — X, whereas the center liner had a depth of about — X. The acoustic resistance of 

8 6 

the center liner was chosen to be twice the resistance of the end liners. The length of the 
center section was made equal to the combined lengths of the end sections. The total duct 
length was chosen to be equal to its diameter. 

The radiation patterns for a point source in this duct are shown in figure 8 and are 
compared with the radiation from a hard- walled duct and from a duct with a uniform liner. 
Radiation patterns for both the flanged and unflanged termination theories are shown for 
comparison. The unflanged duct radiation patterns tend to be 5 to 6 dB lower at 90° than 
the flanged patterns and give results for angles greater than 90°. The admittance for the 
uniform liner is the optimum admittance for a plane wave source in an infinite duct. This 
admittance was given by Rice in reference 10. The results in figure 8 show that this uni- 
form liner is not effective for the point source in a short duct. The three- section liner 
configuration shows substantial reductions in the radiated noise. This single example 
demonstrates the potential usefulness of multisection liners. A complete analytical and 
experimental study of this concept is desirable but is beyond the scope of the present 
paper. 


UNIFORM MULTICHANNELED DUCT (FIGS. 2 AND 9) 

General Equations 

When the duct is divided into many channels, such as the example duct in figure 9, it 
is necessary to write a more general set of transmission and reflection equations than 
equations (17) and (18). For example, right-moving waves at station 6 (fig. 9) may be 
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Figure 8 .- Radiation from source in short duct, f = 1000 Hz; 
c a = 338 m/sec; D = 0,39 m. 



Figure 9 .- Uniform multi channeled duct. 


transmitted to stations 8 and 7 and reflected into planes 6, 5, and 4 as backward- moving 
waves. Similar possibilities exist for other regions which are bounded by planes within 
the duct. Communication between the various terminal planes is also possible as shown 
in figure 9 where right-moving waves at plane 18 are reflected into planes 17, 18, and 19 
and transmitted into planes 1,2, and 3. In general, a matrix may be written for the duct 
mode amplitudes in the form 
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Equation (70) contains all possible interactions of the waves in different planes. Equa- 
tions (17) and (18) for the single channel duct are special cases of the rows of equation (70). 

Axisymmetric Duct 

Splitters and bifurcations .- When an axisymmetric duct is divided into multiple 
channels with circumferential splitters, reflection and transmission matrices must be 
derived which are more general than the ones for the admittance discontinuity; however, 
the general method of derivation, with conservation conditions, is still applicable. 

Consider a duct which is divided by 
splitters as shown in figure 10. At some 
plane in the duct, one set of splitters ends 
and another begins. Again, planes spaced 
a small distance e apart about the dis- 
continuity are used to derive the govern- 
ing equations. Each plane which cuts 
through the splitters defines a number of 
annular regions in which there are modes 

^m ^ r )’ but these modes have no meaning 
except in the range (a) = r = b^. There- 
fore the modes are defined to be zero out- Figure io.- Duct with splitters, 

side of their range, that is 
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xb^ (r) = 0 


if r < a J or if r > 


bi) 


(71) 


The conditions of continuity of axial velocity and pressure (eqs. (42) and (43)) are still 
valid if 



where i = ±j or ±k. Also, since the matching equations (42) and (43) are valid with this 
generalized notation, equations (46) to (53) are also valid. 
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Matrix partitioning . - When there are duct splitters, the mode integral matrix given 
by equation (46) must be partitioned for computational purposes. The partitioning is done 
so that the mode integral submatrices are compatible with the modal subvectors 

^m( X mV r ) in ec l uation ( 72 )- The Partitioning serves to define limits on the integrals 

which generate the submatrices. Since the range of the integral in equation (46) is 
(a = r '= b) and the modal functions associated with any two planes i n and K n are non- 
zero only in the ranges (a^ ^ r § b n ) and (a n ^ r 2 b n ), respectively, the range of the 

integral for the submatrix is^ 11 ^ r ^ b n ) n (a* n ^ r 2 b n )). When this intersection is 

the null set, the submatrix is null. As examples, the modal submatrix integrals for 

3 1 

planes and kg in figure 10 are taken from a to b and the modal submatrices 

for planes jj and kg are null since planes jj and kg do not have a common range 
in r. 

Variable degrees of freedom .- It has been tacitly assumed that the number of 
degrees of freedom is the same at adjacent planes such as j and k in figure 6. This is 
not a necessary requirement for a matching of the acoustic fields and it may be seen that 
equations (47) and (51) are valid even if the orders of the modal vectors at zJ are dif- 
ferent from those at z k . In this case, however, it is not possible to obtain equations 
like equations (49), (50), (52), and (53). Consequently, equations (17) and (18) should be 
replaced by the more general equations (47) and (51) which will permit variable degrees 
of freedom to be used in computational work. 

CONCLUDING REMARKS 

A set of acoustic equations has been derived for the sound field in and radiated from 
multisection ducts. The equations are in matrix form and are based on the use of acoustic 
modes as generalized coordinates. The results of this paper are directly applicable to 
duct systems which are made from a number of joined circular and annular sections. 

Duct systems of this type may be used as a mathematical model of an aircraft engine to 
study the interactive effects of internal noise sources, transmission through complex duct 
liner and splitter configurations, and radiation to the far field. 

There are several advantages to this formulation of the duct acoustics problem. 

The matrix notation used is convenient for programing and digital computation. The 
notation may also be used in variable-area ducts with flow gradients so that the extension 
of the present work to this more general case be done in a routine manner. In addition, 
the matrix formulation may be extended to include ducts with radial splitters. This 
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extension could utilize the acoustic modes in annular sectors and the matrix equations 
would have terms coupling the circumferential modes. 

The matrix formulation of the acoustic equations gives a format for integrating the 
results of several areas of research. Turbomachinery noise sources can be described 
as sources of acoustic modes in an infinite duct. The matrix theory then shows how to 
incorporate these source effects in finite ducts. Transmission through a duct section of 
any type, such as a variable duct or a constant area duct, may be studied without concern 
for the ducts to which it is to be joined since the matrix notation facilitates the joining of 
different duct sections. A variety of mathematical methods may be used: numerical 
methods for very irregular duct sections, perturbation methods for ducts with slowly- 
varying properties, and analytical methods for uniform ducts. The theory of radiation 
from ducts may also be developed using semi-infinite duct models. The matrix equations 
then combine the radiation effects into the finite ducts. Finally, the equations for joining 
different duct sections are identical in form to the equations for transmission through a 
duct section; therefore, all source, duct transmission, and radiation effects are combined 
in a unified format that makes the present theory suitable for programing and computation. 

It has been shown by Lansing and Zorumski (Journal of Sound and Vibration, 

March 8, 1973) that multisection liners can be more effective in reducing radiated noise 
than uniform liners. The present theory shows the two important physical effects by 
which multisection liners reduce radiated noise. These are the transmission effects and 
the reflection effects. In the past, all practical efforts to reduce aircraft noise using 
acoustic treatment have utilized transmission effects. No effort has been made to take 
advantage of reflection effects; however, it has been shown here that they can be impor- 
tant. In a uniform circular duct without flow, the mode reflection matrices at the inter- 
face between two liner sections are shown to be proportional to the differences in the 
acoustic admittances of the liners. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., March 22, 1974. 
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